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We introduce a new Symanzik improved action by adding a 2 x 2 plaquette in sucli a way that the Feynman rules 
in the covariant gauge simplify. We call this the square Symanzik action. Some comparisons with the continuum 
and the standard Wilson action are made in intermediate volumes, where mass ratios are accurately known and 
the precise amount of improvement can be determined. Ratios of the Lambda parameters will be presented, as 
well as partial results for the one-loop improvement coefficients. We discuss some of the intricacies that arise 
because of violations of unitarity at the scale of the cutoff. In particular we show how a field redefinition in the 
zero-momentum effective action allows one to remove scaling violations linear in the lattice spacing. 



1. INTRODUCTION 

We consider here the Symanzik improvement 
scheme Q, which is designed to remove lattice 
artefacts by adding irrelevant operators to the 
lattice action, whose coefficients are tuned by re- 
quiring spectral quantities to be improved to the 
relevant order (on-shell improvement Per- 
turbative calculations, although difficult, are still 
manageable. For Symanzik improvement to work 
it seemed that unreasonably small values of the 
bare coupling constant were required. 

Mean field inspired Symanzik improvement 
H was introduced to beat the bad convergence 
of perturbation expansions in the bare coupling 
constant. In particular the Parisi mean field cou- 
pling 1^ defined in terms of the plaquette expec- 
tation value is seen to improve considerably the 
approach to asymptotic scaling. Despite some at- 
tempts Q no good theoretical understanding for 
this is available. In addition the prescription is 
argued to include tadpole corrections to the coef- 
ficients in the Symanzik improved action, which 
can be seen as a mean field renormalization of the 
link variables on the lattice. Only phenomeno- 
logical arguments have been provided to support 
this. 

One difficulty in testing improvement is how 



to determine to which extent improvement has 
actually been achieved. For pure gauge theories 
standard tests involve restoration of rotational in- 
variance in the heavy quark potential It be- 
comes more problematic when one has to base the 
judgement on carefully extrapolated Wilson data. 

These problems inspired us to consider testing 
improvement for the pure gauge glueball spec- 
trum in intermediate volumes, particularly em- 
phasizing the need to test improvement of scaling. 
In spectroscopy asymptotic scaling is not such an 
important issue since one has to set the scale by 
fixing one of the masses anyhow. The main rea- 
son for considering intermediate volumes (up to 
0.75 fermi across) is that this volume range can 
be accurately described in terms of an effective 
zero-momentum model, nevertheless incorporat- 
ing important non-perturbative features that con- 
tribute to energy of electric ffux. For SU(2) re- 
sults are known both for the continuum limit and 
for the Wilson lattice action, from which precise 
statements on the scaling violations for the mass 
ratios can be made. 

2. SQUARE SYMANZIK ACTION 

As usual, one connects a continuum configura- 
tion with one on the lattice by parallel transport 
of the vector potential along the links 



* Based on the talks "Testing Improvement" and "Hamil- 
tonian from Improved Action" by the last two authors at 
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This allows one to extract the irrelevant higher 
order operators that need to be cancelled in a lat- 
tice action. We introduce the new class of actions 
by adding a 2 x 2 plaquette to the ones considered 
by Liischcr and Weisz, 
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— (co + 20ci- 4c2 + 4c3 + 64c4)^ Tr (P^F^. {x)f 



ca) ^ Tr(I?^F^A(2;)I?,F„A(a:)) + 
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^ ^ Tr((P,F.,)^)+0(a«). 

Sometimes in the literature C2 and C3 are inter- 
changed 0,1 . Here we followed the convention of 
ref. y] and we have taken the liberty of assign- 
ing the coefficient C4 to the 2x2 plaquette. The 
<> imply summing ^ ^ v (7^ A), labelling the 
edges of the plaquette, with the point x attached 
to (say) the lower left corner Q. At tree- level 
only the planar loops are considered, C2 and C3 
acquire non-zero values only at one-loop order, 
but as was shown by Liischer and Weisz ||] field 
redefinitions allow one to put C3 = 0. For the LW 
Symanzik action (c4 = 0) one has cq = 5/3 and 
ci = —1/12 at tree- level, but this does not allow 
for a "covariant" gauge condition that will make 
the gauge field propagator diagonal in the space- 
time indices. The 2x2 plaquette allows one to 
"complete a square" when choosing C4 • cq = c^, 
leading to the gauge fixing functional (z = ci/cq) 
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^z{2 + dl){2+d^))q^{x). (3) 



We decided for this reason to call it the square 
Symanzik action. Here 9^ denotes the lattice dif- 
ference operator d^ip(x) = (p{x + /i) — ^{x). As 
a bonus we note that the condition C4C0 = is 
invariant under multiplicative link renormaliza- 
tion, as they appear in the tadpole improvement 



scheme, allowing one to easily include such factors 
in a perturbative calculation. At tree-level one 
now finds cq = 16/9, ci = —1/9 and C4 — 1/144. 
We have verified that this action satisfies the pos- 
itivity bound Q . It is amusing to see the expres- 
sion for the term in the expansion of the action 
simplify to 
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(behaving as the Symanzik action for cooling |7j). 

At tree-level a tadpole parameter uq modifies 
z = ci/co ~ —1/16 to z = — 1/16wq, whereas 
Co = l/(l-f 4z)^. One finds in the covariant gauge 
the ghost (P) and vector (P^^) propagators to be 

Pik) = ^ 



coEa (4sin^(fcA/2) +4zsin" fc^ 



^(H-4zcos2(fc^/2)) 



(5) 



3. EFFECTIVE ACTION 



One can now perform a background field cal- 
culation to determine the one-loop effective ac- 
tion for the zero-momentum gauge fields. For the 
Wilson action and the continuum this was done 
previously ^ and shown to lead to rather accu- 
rate results. Like for the Wilson case one writes 
U^{x) = e«'^(^)e'=^(*)/^, with N the number of 
lattice sites in the spatial direction (taking the 
number of sites in the time direction infinite) and 
Qfi^x) the quantum field, restricted to non-zero 
(spatial) momentum, to be integrated out. This 
choice on splitting off the quantum component of 
the lattice field yields a particularly simple back- 
ground gauge fixing function 



9f 



/Co 



(7) 



"covariantizing" the difference operator to 
DM^) EE e"^(*)/^^(a; + /i)e-"^W/^ - ip{x). 

3.1. The SU(2) effective potential 

The one-loop calculation greatly simplifies for 
an abelian constant background field as this al- 
lows one to diagonalizc the propagator with re- 
spect to the isospin neutral and charged decompo- 
sition of the gauge and ghost fields. The momenta 
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in the background field c = ^iCa^ are shifted, 
k-^k + sC /N , where s = for the neutral isospin 
component and s = ±1 for the two charged com- 
ponents. It is not hard to find the eigenvalues of 
the fluctuation operators for the ghost and gauge 
fields 

\gh{k) = VQ^^4sin2(fc,/2)(l+4zcos2(fc,/2)), 

Ap(fc) = V^(l + 4zcos2(A:^/2))Agft(fc). (8) 

These eigenvalues can be written as products of 
factors 4sin^(fco/2) -t-w^, where the uja can occur 
in complex conjugate pairs at spatial momenta 
close to the edge of the Brillouin zone. As it is 
well-known that the sum over fco for one such a 
factor can be performed explicitly, it is not sur- 
prising we can find a closed expression for the ef- 
fective potential, as a sum over the appropriately 
weighted logarithm of the eigenvalues 



Vr{C)^N}^\^^\og (H-4zcos^ 



27:711 + Ci 



2N 
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with tj2 = 4^.sin2(fc^/2) (l + 4zcos2(fc,/2)), and 
k= {2Trn + C)/N. This effective potential, nor- 
malized to V(0) = 0, is plotted for uq = 1 in 
fig. 1 as compared to the result for the Wilson 
action (z = 0) and for the continuum (A^— >oo). 
Although this effective potential is not spectral, 
since near C G 27^71? the adiabatic approximation 
for integrating out the "charged" zero-momentum 
modes breaks down, one sees that improvement is 
quite efficient in removing scaling violations (only 
scaling violations to fourth order in the lattice 
spacing a = 1/N remain). At A^ = 6 we can 
not distinguish the result from the continuum at 
the scale of this figure. One might even fear that 
choosing uq 7^ 1 makes the agreement worse. 

3.2. The Lambda ratios 

One can proceed as in the Wilson case with 
computing the one-loop coefficients for the effec- 
tive action In particular this provides for 
N 00 the renormalization of the bare lattice 




Figure 1. The effective potential for a constant 
Abelian background field Ci = ^iCa^. The full 
line represents the continuum result. The lower 
two dashed curves are for the square Symanzik 
action with A" = 3 and 4. The upper three dotted 
curves are for the Wilson action with A^ = 3, 4 
and 6. 



coupling at fixed physical volume. As Lorcntz in- 
variance is replaced by cubic invariance (both by 
the lattice discretization and the periodic bound- 
ary conditions), two independent determinations 
of the Lambda ratios can be extracted from the 
effective Lagrangian 



i(l- 

2 \ 2 
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dt 



+ i{^ + a2){Ft^\v,{c). (10) 



Here g'"^ = g^'^ - ll\Qg{N) /12t:'^ is kept fixed 
while taking the continuum limit. In this limit 
ai and a2 differ by identical finite amounts from 
what was found for the continuum and the Wil- 
son action. This finite difference allows one to ac- 
curately compute the Lambda parameter ratios. 
These ratios were also computed using the heavy 
quark potential method of ref. [8], which allows 
one in addition to obtain the result for SU(3). 
For the ratio of the square Symanzik action to 
the Wilson action we find. 



As2/A,, 



4.0919901(1) forSU(2)\ , . 
5.2089503(1) for SU(3) ) ' ^ ' 



3.3. One-loop improvement 

Tests of tadpole corrections to variant tree- 
level improved actions have been performed be- 
fore |10|. However, we remind the reader that 
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there is only one pure gauge improved lattice ac- 
tion that was computed to one-loop order ||]. 
It is our aim to bring the square Symanzik ac- 
tion to this same level. In principle this pro- 
vides a way of testing to what extent the suc- 
cess of the tadpole improvement depends on the 
choice of action. Independently it is a useful 
check on the consistency of the Symanzik im- 
provement scheme with its inherent redundancy 
in choosing the lattice action to cancel scaling vi- 
olations. For these perturbative calculations of 
the one-loop corrections c'^ one follows the well 
established route of using the twisted finite vol- 
ume spectroscopy As a normalization con- 
dition on the definition of the coupling constant 
one imposes Cg -t- 8c']^ -I- -I- I6C4 = 0. Requiring 
the physical mass of the lowest state to have no 
quadratic scaling violations to one-loop order for 
the square Symanzik action leads to 

, _ ,,.,_( -0.00838(1) for SU(2) 
ci C2 + 4C4 - 1^ -0.01545(2) for SU(3) 

As an independent check this combination was 
also extracted (at higher accuracy) from the 
heavy quark potential. In addition the on-shcU 
three point coupling extracted in the twisted fi- 
nite volume allows one to find a value for 
36(c'j — c'2 + 4C4) -I- 8C2. This computation is 
rather involved and still in progress. Note that 
C4 appears in the combinations c'l -\- ic'^ and 
Cq — I6C4, as is also dictated by eq. (2). There- 
fore as was to be expected C4 is a free parame- 
ter. It need not, but can, be fixed by requiring 
C4C0 = c\ to one-loop order. Finally we quote the 
result for the single plaquette expectation value: 
= 1-0.35878 •g§(A^-iV-i)/4. 

4. MONTE CARLO RESULTS 

In the intermediate volume context we can con- 
sider at most lattice spacings up to 0.25 fermi 
(asking to be absolved for pushing in this direc- 
tion), as the larger plaquettes that appear in the 
improved actions require the volume to be at least 
three lattice spacings in each direction. Despite 
the appearance of unitarity violations at the scale 
of the cutoff, due to the non-local nature of an im- 
proved action, the intermediate volume physical 
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Figure 2. SU(2) Monte Carlo data for mass ratios 
in a small volume on a lattice of size 4^ x 128, us- 
ing the Wilson action (crosses from existing data 
of Michael), the LW Symanzik action and our new 
square Symanzik action. The lines give the an- 
alytic results: full for the continuum and dotted 
for the standard lattice action {N = 4). 



masses remain small enough in lattice units to ex- 
tract them from the decay of correlation functions 
in the time direction in the usual way. In larger 
volumes and coarser lattices the latter problem 
was dealt with by taking the lattice much finer in 
the time direction pl] | , using the asymmetric cou- 
plings well-known from finite temperature stud- 
ies [Q. Here we will consider only Monte Carlo 
data at 0.018 fermi {P = 4/gg = 3 for the Wilson 
action, for which we compare our data to those 
by Michael |Q). Odd as it may seem, this is 
where the scaling violations for a lattice of 4 lat- 
tice spacings in the spatial directions are largest 
within the finite volume spectroscopy The 
data corresponding to the LW Symanzik action 
is represented by the triangles and for our new 
square Symanzik action by the squares. In both 
cases we used tree-level improvement only. The 
improvement is significant. For the LW Symanzik 
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action the data is within two sigma of the contin- 
uum values. The results seem to indicate that 
the square Symanzik action is somewhat less ef- 
fective, although the difference is not significant. 
A comparison at coarser lattices will be more in- 
teresting as one should like to see, as advocated, 
tadpole corrections to further improve the results. 
For this purpose we present data elsewhere at a 
lattice spacing of 0.12 fermi. 

5. HAMILTONIAN 

5.1. Toy model 

A well-known problem of improved actions is 
that the transfer matrix is not hermitian Jlif . 
This is easily seen to be related to the next- 
to-nearest neighbor couplings in the time direc- 
tion. We will illustrate things here by means of 
a simple one dimensional problem. For the ac- 
tion S{x) = J2tix(t+^)-xit))^ / {2g^a)+aV{x{t)) 
the partition function at finite Euclidean time 
(T = aN) can be exactly rewritten in operator 
form 



Vxe 



-S{x) 



Tr(e 



where K — —^g^{d/dx)^ . The Hamiltonian read 
off from this equation is only determined up to a 
unitary transformation. To lowest order one finds 



H = K + V 



-[V,[K,V]]+0{a^). 



(14) 



Note that [V, [K, V]] = g^V'{xf. 

Next improve the kinetic term {x{t+l)— x{t))'^ 
by 4(a;(i-Hl)-a;(i)) 2/3 - {x{t+2)-x{t)f 112. One 
finds that the propagator factorizes as P{k) = 
(P_ (fc)-P+(A:))/Z, where P± \/c) = 4 sin2(ifc>f w^, 
Z = ^l-aW/a, uj'i = 6(1 ± Z) and = 
g^V"{Q). This explicitly exhibits the unphysical 
pole mentioned before with masses 12 /a 

at the scale of the cutoff. They are not harmful 
for low-energy behavior [p^ . It would perhaps be 
misleading to associate the spurious poles with 
ghosts as they do not just occur in loops. Ver- 
tices do not preserve ghost number. Nevertheless 
we expect their contribution to low-lying states 
to be suppressed in a way similar to the influence 
of virtual processes due to heavy particles. 



Let us introduce the following field redefinition 
best expressed in the Fourier representation 

a^g^ dVix) 



:{k)^x{k)Jl + J^k^- 



24 dx{-k)' 



(15) 



where as usual fc = 2 sin(fc/2). When substituting 
this non-local transformation in the action we find 



k 

= E 



k^{l + J,k^)\x{k)\^ 



2ag^ 



-aVix) 



(16) 



e\xik)i 

2a.g2 



+ aV{x) + ^r{xr+0{a'), 



for which it is assumed that k ~ 0{a). We 
note that after the field redefinition, ignoring the 
0(0^) corrections, the action is local in time and 
one obtains H = K+V+0{a'^) from eq. (14). 

However, interactions will give rise to a non- 
trivial Jacobian under this change of variables, 
J{x) =detidx{k)/dx{k')), or 

..d'Vix) 



(17) 



J'(^)=dct(^l + ^5t9o-^aV^ 

where we took the liberty of modifying the 0{a*) 
terms in the operator whose determinant is to be 
evaluated. We could likewise define the transfor- 
mation such that the Jacobian is given as above, 
although this leads in multi-dimensional cases to 
non-integrable transformations. 

Remarkably one can rewrite this Jacobian up 
to an X independent factor as 



^^(-)=det(l-aVP+[^ 



(18) 



Its contribution to the partition function Z can be 
interpreted as the effective action in a background 
field calculation, with the propagator truncated 
to the unphysical branch, albeit to lowest non- 
trivial order in the lattice spacing. One easily 
verifies that for V{x) = Aa;^ this Jacobian gives 
rise to a mass correction linear in the lattice spac- 
ing, which was initially discovered by computing 
the mass gap to first order in A from the Feyn- 
man rules for the improved action of this simple 
model. As the model is quite similar to the ef- 
fective action we discussed before, we were quite 
puzzled by this result and it prompted the above 
derivation. 
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5.2. Gauge model 

Indeed, taking the results of sect. 3.1 we can 
compute easily part of the effective action for the 
zero-momentum gauge fields. To obtain the effec- 
tive potential that is valid near C — 0, where the 
tree-level potential is quartic in the gauge fields, 
one restricts the sum to n ^ and replaces Ci 
by Ti = -^/^2Trcf in eq. (9). The result is de- 
noted by Vi{r). An accurate description of the 
full effective potential to one-loop order is given 
by 

Vi (c) = Vi ir)+a3rfF^^^^+airlFtj Vagdet^c (19) 

As can be extracted from the zero-momentum 
part of eq. (9), 7i(Af) =71(00) + A(V3-l)/iV-|- 
0{1/N^), where 71 is the coefficient of r ^ in the 
effective potential. Generalizing the analysis of 
the toy model to the situation at hand one finds 
from the Jacobian Sjji — —^^/S/N. The miss- 
ing piece is provided by the non-triviality of the 
Haar measure for integration over the background 
link variables 

ShVi{c) = -2NY,^og[2Nsm{n/2N)/n], (20) 

i 

Indeed one finds Sh'Ji — t^A^^^- Both the Jaco- 
bian and measure contributions compensate for 
the scaling violations linear in the lattice spacing 
and with it Vi(r) becomes free of scaling viola- 
tions to third order in the lattice spacing. 

Furthermore, rescahng c with (1 + J^g'^^i/N'^) 
removes to a high degree of accuracy unwanted 
scaling violations in ai and a2- More sur- 
prising was to find that the field redefinition 
5c^ = -2g'^\og{N)VlF^i/{2^nNf is required to 
remove \og{N)/N'^ scaling violations in and 
as. As a non-trivial check the one- loop coeffi- 
cient ao, in front of the term ^(Sq^oc")^, was 
computed. Its log(A^)/A^^ term combines after 
the above field redefinition with the tree-level co- 
efficient of l/(12(7gA^^) in precisely the right way 
to renormalize the coupling constant. 

Remaining 0{N~'^) scaling violations will and 
can be cancelled by the one-loop improvement co- 
efficients. Due to the unfortunate mixing with co- 
efficients that are not easily accessible in lattice 
perturbation theory |^], we cannot at present get 
at these one-loop improvement coefficients along 



this route. As we have seen, using an effective 
action one imposes improvement only up to field 
redefinitions (or up to unitary transformations at 
the Hamiltonian level) . This is more difficult than 
computing a few spectral quantities, but has the 
obvious benefit of manifestly improving infinitely 
many levels at the same time and would be a very 
non-trivial check on improvement indeed. 
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